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Abstract 

A Grobner basis-based algorithm for solving the Frobenius Instance Problem is presented, and 
this leads to an algorithm for solving the Frobenius Problem that can handle numbers with 
thousands of digits. Connections to irreducible decompositions and Hilbert functions are also 
presented. 



1. Introduction 

Let pi, . . . ,p n be relatively prime positive integers and let p = (pi, . . . ,p n ). An in- 
teger is p-representable if it can be written as v ■ p for some v 6 N". Determining 
p-representability is known as the Frobenius Instance Problem, and we present a new 
Grobner basis-based algorithm that solves it PI Our algorith m differs from the c l assica l 
way to s olve int e ger p rograms using Grobner bases due to Conti and Traversol (1991) 
(see also Pottier ( 1994 )) by not adding any auxiliary variables to the problem. 

The Frobenius Number f* is the largest integer that is not p-representable, and such 
an integer exists by Proposition 2 below. E.g. if p = (6, 10, 15) then /* = 29. 

The Frobenius Problem i s to compute the Frobenius Number f*. A recent algorithm 
due to lEinstein et alj (|2007l ) can solve Frobenius problems even if the pi have thousands 
of decimal digits. Here we describe a simpler variant of that algorithm that performs 
better. 

The algorithm f irst computes a Grob ner basis and then determines the Frobenius num- 
ber based on that. lEinstein et al.l (|2007l ) use a novel algorithm based on the Fundamental 
Domain that in effect computes this Grobner basis, and they report that this is much 
faster than using Buchberger's algorithm in their implementation. 



However, the benchmarks in section 7 show that the program 4ti2 [4ti2 team (2006)], 
which implements Buchberger's algorithm in a special case, outperforms the Fundamental 
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Domain-based implementation from lEinstein et al.l (|2007t ). The record for random^ with 
11 decimal digits was n = 11, but we can now reach n = 13. Performance is also improved 
on smaller examples. 



We show that the second step of the algorithm of lEinstein et al.l ( 20071 ) can be rephrased 
as the computation of the irreducible decomposition of a certain monomial ideal. The 
final step of the algorithm is to maximize a linear function over the decomposition, and 
we show that this essentially computes the index of regularity of the ideal. We define 
these terms when they are needed. 

See iJensen et al.l (|2007l ) for a geometric version of some of the ideas in this paper in 
terms of maximal lattice free bodies. That paper is joint work with Niels Lauritzen and 

Anders Nedergaard Jensen. 

We refer to the book by Ramirez Alfons in ( 2005) for more background on the Frobe- 
nius Problem. We wish to thank Anders Nedergaard Jensen, Niels Lauritzen, Daniel 
Lichtblau and Stan Wagon for helpful discussions about Frobenius numbers. 



2. Preliminaries 

Let ei G N n , i = 1, . . . , n, be the vector whose entries are all zero except that there is 
a 1 at position i. The p- degree of a vector v € Z n is v ■ p. If v € Z n then define v + G N™ 
as below and define v~ := (— v) + . 

Vi, for Vi > 
0, for Vi < 

We will need to consider the lattice ideal Ic defined by 



x v+ - x l 



v G Z n and v ■ p = 



We will compute a Grobner basis G of Ic- The term order < we will use first considers 
the p-degree of the exponent vector of a term and then the reverse lexicographic order 
where x\ < X2 < ■ ■ ■ < x n . 



Example 1 I f y = (2,3) then I = x^ < x^ < x^ < x {0 ^ < x^ 



4^5 



We refer to Cox et al. I (|l997h for more details. 
Proposition 2 Only finitely many integers i £ N are not p-representable. 
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Proof. As pi, . . .,p n are relatively prime, iterated use of the Extended Euclidean Al- 
gorithm provides us with a vector v G Z™ such that v ■ p = 1. Let m := min™ =1 Vi and 
define u :— v + pi|m|(l, . . . , 1). Then u + iv G N n for i = 0, . . . ,p% — 1 and therefore 
(u + iv) ■ p = u ■ p + i are p\ consecutive p-representable numbers. □ 



3. An Algorithm That Solves The Frobenius Instance Problem 

We claim that the following algorithm determines if t G N is p-representable. 
Step 1: Compute an a G Z™ such that a ■ p — t, a\ < and Oj > for i = 2, . . . , n. 
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Step 2: Divide x a — x a by G giving remainder x w (x c — x c ) for some w G N™ where 
x c < x c . 

Step 3: Then t is p-representable if and only if c G N™. 

Step 1 We first find an a G Z™ such that a ■ p = t. One way to do this is to use the 
Extended Euclidean Algorithm iteratively to find a b G Z™ such that 6 • p = 1 and 
then let a := t&. As (— 5Z™ =2 P«'Pi' ■ ■ ■ 'P 1 ) nas P-degree zero and the sign pattern 
(— , +,•■•, +), we can assume that a also has this sign pattern by adding a sufficiently 
large multiple of this vector to a. 

Step 2 This step requires knowing the Grobner basis G, and computing G is the most 
time consuming part of the algorithm. Once G has been computed the polynomial 
division itself is comparatively fast. 

Step 3 Observe that the division algorithm ensures that c will have no more negative 
entries than a does, so c is negative at most in the first coordinate. As x c is not 
reducible by G, the following lemma tells us that if t is p-representable then c\ > so 
that c is a p-representation of t since c ■ p — a ■ p = t. 

Lemma 3 Let a G Z" such that a.j > for i — 2, . . . ,n and a\ < 0. If a ■ p is p- 

representable then there exists a g G G such that in<{g)\x a . 

Proof. As a ■ p is p-representable there exists a b G N™ such that a ■ p — b ■ p. Letting 
d := a — b this implies that d ■ p = whereby h := x d — x d G Ic- Thus there exists a 
g G G such that in<(g)|in<(/i). We will prove that in<(g)|in<(/i) = x d \x a . 
in<(h) = x d+ : d + and d~ have the same p-degree and d\ = a\ — b\ < 0. 

x d+ | x a+ . This follows f rom beW 1 . □ 

Note that < can be replaced with any term order that first considers the p-degree of 
the exponent vector and then the reverse lexicographic order on the first variable. 

4. An Algorithm That Solves The Frobenius Problem 

The general idea of the algorithm is that we can represent /* by a certain vector that 
has p-degree /* , and that this vector has certain properties (see Proposition 4) . It turns 
out that only finitely many vectors have these properties, so we can look through all of 
them, and then the one with maximal p-degree will be the vector that represents /*, i.e. 
it will have p-degree equal to /* (see Proposition 5). 

Proposition 4 spells out the properties mentioned above. 
Proposition 4 Let c G Z™ be the vector resulting from running the algorithm from 
Section 3 on t = f* . Then the following holds. 
(Ml) c\ = —I and a > for i = 2, . . . , n. 
(M2) x c+ cannot be reduced by any element ofG. 
(MS) x' c+e ^ + can be reduced by some g t G G for i = 2,...,n. 

Proof. (Ml): It holds by construction that c\ < and that Cj > for i = 2, . . . , n. 
Let c' := c + e\. Then x^ c ) + = x c+ is not reducible by G and d ■ p is p-representable as 
d ■ p > /*. Then Lemma 3 implies that d x > whereby c\ = — 1. 
(M2): This holds by construction. 



3 



(M3): We see that (c + e.;) • p is p-representable as it is strictly larger than /*. Thus 
Lemma 3 provides a ft G G such that in< (<7i ) |x^ c+ei - )+ . □ 

Let M p be the set of vectors c E Z n that have the properties (Ml), (M2) and (M3) 
from Proposition 4. The idea is to compute M p and then use Proposition 5 below to find 

/;• 

Proposition 5 The following holds, 
(i) M p is finite. 

(ii) If a G M p then a ■ p is not p-representable. 
(Hi) f* = max {a ■ p\a G M p } 

Proof, (i): Let a 6 M p and i 6 {2,...,n}. Let g e G such that m<(g)\x ( - a+e '^ + . As 
in< (g) does not divide x a+ we can infer that + 1 is the exponent of Xi in in< (g) . Thus 
there are at most \G\ possibilities for what ai can be. 

(ii) : Lemma 3 shows that this follows from (Ml) and (M2). 

(iii) : Proposition 4 shows that there is an a £ M p such that a- p — f*. This and part 
(ii) above shows what we need. □ 



We can compute M p as follows. We saw in the proof of Proposition 5 that if i G 
{2, . . . , n} and a G M p then + 1 is the exponent of Xi in in<(ft) for some gi G G. Thus 
we can run through all possible values of for i = 2, . . . , n and only keep those a that 
have properties (Ml), (M2) and (M3). 

This algorithm is easy to understand and implement, but it re quires us to look throu gh 
up to \G\ n 1 possibilities. The External Corner Algorithm from lEinstein et ah ( 2007 ) is 
a more efficient algorithm for computing M p which usually dramatically reduces the 
number of possibilities that need to be examined. 

In Section 5 we will define the irreducible decomposition of a monomial ideal, and we 
will prove t hat M r, corresponds to the irreducible decomposition of the initial ideal of 



Ic- Round ( 2007 ) shows how an algorithm that is very similar to the External Corner 



Algorithm can compute irreducible decompositions of monomial ideals in general in much 
less time than the best competing programs. 



5. A Connection To Monomial Irreducible Decompositions 

We need a few more definitions. Let / be a monomial ideal and define the function 
4> by (f> (v) = (x^' \vi > 0) for v G N" except that <f> (1) = (1). An ideal of the form <p (v) 
is called irreducible and the irredundant irreducible decomposition of / is the unique 
minimal subset D G N™ such that / = f] veD <fi (v). Thus the irredundant irreducible 
decomposition of I := (a;f , X1X2) is {(1, 0), (2, 1)} as / = (xi) fl {x\, x 2 ). 

An ideal is artinian if there exists a t G N such that x\ G / for i = 1, . . . , n. Note 
that in<(/c) is artinian if we first project out the variable xi, since x? 1 — x\ l G Ic and 
therefore in^x^ 1 — x\*) = x\ x G in<(/c) for i = 2, . . . , n. 

We claimed in Section 4 that M p corresponds to the irreducible decomposition of the 
initial ideal of Ic, and Proposition 6 proves this claim. 

Proposition 6 The set M' :— {(ai+1, . . . , a n +l)\a G M p } is the irreducible irredundant 
decomposition D of in<{Ic)- 
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Proof. Let a G Z" and let a' := a — Yl7=i e i- Consider the following statements. 

(1) aeD 

(2) di = 0, <n > 1, m<(Ic) C (a) and in<(7/;) $2 ( a + e i) for i = 2, . . . , n. 

(3) < = -1, a' t > 0, x( Q ') + £ in<(7 £ ) and a;( Q '+ e + e i n <(7 £ ) for i = 2, . . . , n. 

(4) a' G M p 

We will prove that these statements are equivalent. 

(1) <^> (2): The initial ideal of Ic contains no monomial that is divisible by x\, and 
this implies that a\ = for a G D. Now that we have handled the first entry, we can 
project out x\ and thereby get an artinian ideal. 

When working with artinian ideals the elements of the decomposition consists of vec- 
tors without zero entries, and if a has no zero entries, then (f>(a + a) C <p (a). 

To get the irredundant irreducible decomposition of an ideal we write the ideal as an 
intersection of irreducible ideals that are as small as possible, and this is exactly what 
(2) expresses since the projected ideal is artinian. 

(2) <^> (3): By Lemma 7 below. 

(3) <^> (4): By definition. □ 

Lemma 7 Let m G N" and let I be a non-zero monomial ideal. Then x rn G I if and 
only if I %4>(m + YJi=i e i)- 

Proof. Let a e I be a monomial. Then a\x m if and only if a ^ § (m + Ym=i e *)- n 

6. A Connection To The Hilbert Function 

The p-weighted Hilbert function HF/:N — > N of a monomial ideal I is defined such 
that HF/(£) is the number of monomials x m £ I where m has p-degree t. 
Proposition 8 HF^ n< ^^ (t) is equal to 1 if t is p-representable and otherwise. 

Proof. HF m< ^^ (t) < 1: Let mi,m 2 be two vectors of p-degree t. Then x mi — x™ 2 G Ic 

whereby the initial term is in in<(/c). Thus at most one of x mi and x" 12 is not in in<(l£). 
HFj n ,j s (t) = 1 => representability: If HF m< ^ j£ ^ (t) = 1 then there is a monomial 

x m £ in<(/£) where m has p-degree t. Thus t is p-representable. 

representability => HFj n<( - Ic ^ (t) = 1: Run the algorithm from Section 3 on t. The 

resulting vector c G N" is such that x c £ in< (Ic). □ 

We can infer that HF m< ^ j£ ^ (/*) = and that HF m< ^ £ ^ (t) = 1 for all integers 
t > f*. The integer at which the Hilbert function becomes equal to a polynomial is 
known as the index of regularity, so in this case the index of regularity is /* + 1 . 

The elements of M p correspond to the maximal monomials outside of in<(/c) accord- 
ing to divisibility (disregarding the first variable), and what the algorithm from Section 
4 does is to maximize the dot product with p over M p . This amounts to maximizing the 
dot product over the vectors m such that x m+ei £ in<(/£) and mi = —1. 

We can multiply any monomial not in in<(i£) with x\ and get something still not 
in in<(/c), so anything that goes on in the first variable does not prevent the Hilbert 
function from becoming a polynomial. Thus the algorithm can be interpreted as finding 
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the maximal point that prevents the Hilbert function from becoming a polynomial, which 
is to say that it computes the index o f regularity. 



We conclude that the algorithm of lEinstein et al.l (|2007l ) can be interpreted as com- 
puting an index of regularity. 

7. Benchmarks 



Round §006) has written an implementation called Frobby of the algori thm described 



i n thi s paper, and here we compare Frobby to the implementation of lEinstein et al 



(|2007l ). Figure 1 displays the collected data. "Intractable" means intractable according 
to the authors of that software package. Note that the most time consuming part of the 
algorithm usually is to compute the Grobner basis. Frobby and Mathematica are the 
only programs that can handle input nu mbers p; as la rge as those in figure 1 



ly programs that can handle input numbe rs Pi as la rge as those m hgure 1. 
Frobby uses the program fplll due to IStehl to obtain an LLL-reduced lattice 



basis and the n compute s the G robner basis of the lattice ideal Ic from that using the 
program 4ti2 4ti2 teaml ( 20061 )]. Frobby then computes the Frobenius n umber by com - 



puting an irreducible decomposition of in<(i£) using the algorithm due to lRound (|2007l ). 
which is similar to the External Corner Algorithm. Frobby is written in CH — h and the 
source code is available un der the GNU General Public License (GPL) 



The implementation of lEinstein et al.l ( 20071 ) is available as a part of Mathematica 



and is written mostly in C. The major difference from Frobby is that the Grobner basis 
is computed using a different algorithm than Buchberger's. 

All the inputs were r andomly gener ated using genuinely random radioactive decay via 
the service provided by IWalkerl (|2006l ) except the n = 11 input which was provided by 
Stan Wagon who pseudo-randomly generated it using Mathematica. We wish to thank 
Daniel Lichtblau for carrying out the Mathematica benchmarks. 

All the benchmarks were run on machines with a 3.0 GHz Pentium 4 CPU with 1 GB 
RAM except the Mathematica benchmark on the n = 11 input which was run on a 3.2 
GHz Pentium 4. 
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Fig. 1. The benchmark data. 
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